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Abstract 

The generalized Bretherton equation is studied. The classifica- 
tion of the meromorphic traveling wave solutions for this equation is 
presented. All possible exact solutions of the generalized Brethenton 
equation are given. 



1 Introduction 

In recent papers ilii2j, we proved two theorems for possible presentation of 
solutions for nonlinear ordinary differential equations on the complex plane. 
These theorems allow us to have a classification of exact solutions for nonlin- 
ear differential equations and lead to a new method for constructing mero- 
morphic exact solutions of some nonlinear ordinary differential equations. 
In the case one branch for expansion of the general solution in the Laurent 
series we obtain the known methods for finding exact solutions which were 
developed in the last years [SHIlj. However for two or more branches of ex- 
pansions for the general solution in the Laurent series we can look for new 
exact expressions for exact solutions of nonlinear differential equations. 

The essence of approach from papers [TII2] consists in looking exact so- 
lutions of nonlinear ordinary differential equations by comparing the Lau- 
rent series for the general solution of ordinary differential equation and the 
Laurent series for functions that can be a solution to the original equation. 
Taking this approach into account the meromorphic exact solution of third 
order differential equation were found in [T]. We also apply our method for 
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constructing the exact solutions of the Kawahara equation expressed via the 
Weierstrass elhptic function in [2]. 

In this paper we look for all possible traveling wave solutions of nonlinear 
evolution equation which is the Bretherton equation and takes the form 

Uu + aU:rx + Pur.r.r.r. + Su"' + 'ju'' = 0, n>m, u^l (1) 

Using traveling wave u{x,t) = y{z), where 2; = a; — Co t in ([1]) we have 
the nonlinear ordinary differential equation 

(a - C',) + /3 y,,,, + Sym + ^yn = Q (2) 

Integrating equation ([2]) once with respect to z we obtain 

v! + y... - 1 yl) + ^ y""^' + ^ y""^' + 0^ = 0, (3) 

/ 2 m+1 n+1 

where C2 is constant of integration. 

In work [2] the full classes of meromorphic traveling wave solutions of ([2]) 
were have obtained at m = 1 in the cases n = 2, n = 3 and n = 5. 

The aim of this work is to search for meromorphic the traveling wave 
solutions of Eq.([T]) at any n and m = 1. We show that the elliptic solutions 
exist at 77. = 2, n = 3, n = 5 and at n = |. We also study equation ([1]) in the 
case n = 5atm = 2,3,4. To achieve our goal we will use approach suggested 
in [nig. 

The outline of this paper is as follows. In Section 2 we show that Eq.(IT]) 
has elliptic traveling wave solutions at m = 1 if and only if n = 2, n = 3, 
n = 5. In Sections 3 and 4 we obtain meromorphic traveling wave solutions 
of Eq.© at m 7^ 1. 

2 Exact solutions of Eq. (I3j) at any n in the 
case m = 1. 

Consider the generalized Bretherton equation in form ([2]) . Denoting a — Cq = 
a we have 

(ryzz + /3yzzzz + Sy + -fy'' = (4) 
Without loss of the generality we take 

P = -l, 7 = 1 (5) 
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From (jlj) we have 



(y Vzz - Uzzzz + 5 ?/ + = (6) 
Let us prove the following theorem. 

Theoremin the case <t ^ and 6 there are elliptic solutions of 
equation IQ if and only if n = 2, 3, 5. 

To prove this theorem we use the approach by paper |15j , which allows us 
to remove arbitrary exponent n from Eq.([n]). Let us introduce the Backlund 
transformations between solutions of Eq. ([6]) and transformed equation which 
has solution with pole of the first order. 

Taking into account the transformations 

Vz 

w = —, 

y (7) 

from Eq. we have 

3w (n — 5)wl + (Saw + 6 nw'^ — 10 — 10 Wzz — naw) Wz+ 
+6w — n6w + nw^ — naw^ + Anw'^Wzz + nwrn^zz — + a (8) 

+aWzz - 5 WzzzW - 10 WzzW'^ - Wzzzz = 

We suppose that solution of Eq. (jS]) can be presented in the form of Lau- 
rent series in an neighborhood of movable pole z = Zq 



oo 

w[z) = 

k=0 



J2akiz-zo)''~P, p>0 (9) 



The general solution of Eq.Q has four different Laurent expansions in 
a neighborhood of the movable first order [p = 1) pole z = zq. Taking the 
autonomous of equation into account we can take zq = 0. Substituting 
series ([9]) into Eq. (|8]) we have four following values of 

4^) = 1, af = 2, = 3, 4^) = -^. (10) 

The necessary conditions for the existence of the elliptic solutions of equa- 
tion Eq.(|H]) we can find from the following equations [T|I^[TB|[T7] 

4^^+a[,') = 0, = 1,2,3), (11) 

+ = 0, (z,j = 1,2,3), t^j, (12) 
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a« + + + = 0. (13) 

Solving equations f ITT]) . f lT2]) and ffT^ we obtain the following values n for 
existence of the elliptic solutions of Eq.(|H]) 

5 9 7 

n=-, n = -, n = 2, n=-, n = 3, n = 5 (14) 
3 5 3 

Cases n = 2, n = 3 and = 5 were investigated in |j2j. Let us consider 
the cases of = |, n = | and = |. 

3 ' 5 6 

Note that we should consider the Laurent series for solution of Eq. ([H]) 
corresponding to coefficient Oq''^ [15j. 

2* The COr± \::;ojj winj-ixi^ vdiu-^o wi U/q 

= 1, = 2, = 3, 4''^ = -6 (15) 



Consider the case n = |. The corresponding values of following 



We see that elliptic solution may be obtained if and only if we use the solution 
corresponding to combination of the four asymptotic expansions. 

Following to the method [T],i2] we construct the formal Laurent expansion 
of Eq.([H]) solution corresponding to coefficient a^^'' = —6 



6 a 50(5+107or2 

w = z z - 

z 50 615000 



(16) 



129150000 



Let us note that coefficient qiq in the Laurent series ( 1T6|) is arbitrary. 
In accordance with the method PQ[2] the general elliptic solutions has the 
following form 

2(p-A) + 2(p-A,) + 2{p-A,) 

where p{z,g2,g3) is the Weierstrass elliptic function. 

Now we must compare expansion (1161) with the Laurent expansion for 
expression (1171) . The Laurent expansion for (1171) may be obtained by using 
addition formulas for p and pz- The Laurent series and addition formulas for 
p and pz are given in [18]. Also the Laurent expansion for expression ([1] 
may be obtained with help of symbolic computation software Maple. 
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Comparing expansion flTBl) with the Laurent expansion for expressions 
(HZ]) we get 

Ai = ^2 = A3 = — — , Bi = , d 



300' 2 ' 625 ' /ION 
_ 3617^^ ^ ' 

~ 7500' ~ ~ 3375000' ~ 271367971875000000000000000 
We see that there is constraint on the coefficient aig. Otherwise solution 

in the form f[T7|) does not exist. 

However at these values of invariants g2,g3 given in (|T8i) the Weierstrass 

elliptic function is degenerated to trigonometric function 

P ^. -=1 4 tanh^ 1 44) (19) 



7500' 33750007 150 V 2 I 10 _ 

Substituting formulae ffT^ into f lT7|) and taking f lTS]) into account we have 
solution of the Eq. in the from 



w = -^v^tanhj^^l (20) 

We can see that solution (|T7|) is degenerated to the simple periodic solution. 
In this case the exact solution of the Eq.([6]) has the following form 

y = a-^ cosh z (21) 

^ 15625 V 10 / ^ ^ 



(22) 



In the case of n = | values of following 

.(1) _ 1 .(2) _ o ^(3) _ Q ^(4) _ c; 



Therefore elliptic solution may be obtained if and only if when we consider 
solution corresponding to combination of ag^-*, Og^'', . 

Accordingly with method from [li[2] the solution can have the following 
form 

As in the previous case we construct the formal Laurent expansion of 
5 5(T / 137(t2 5 



Eq.([8]) solution at = | corresponding to coefficient a^^"^ = —5 



w = z — I 1 I z'^ 

z 222 V 492840 750 



a (72775 + 323084 5) 5 ^3 

— z + . . . + ai4 2 + 



(24) 



19146834000 



where coefficient ai4 is an arbitrary constant. 

The Laurent expansion for fl2^ can be calculated using formulas from 
textbook [18] or using Maple command. Comparing the expansion fl24l) with 
the Laurent expansion for function (l23l) we have 

3^2 

,3 = -^!^. = A (25) 

5476 ' ^ 1367631' ^ 4107' 

ai4 - 



11083927002941913120 

We obtain constraint on the coefficient 014. Only in this case Laurent series 
corresponding solution (12^ exist. 

However at the values of parameters elhptic solution (12^ degenerate 
to the following simple periodic solution 

w = -S^^tanh (^y^ (z - zo)^ (26) 

In this case the exact solution of the Eq.([n]) has following form 

420^/4 1369^/4 . ,„ , 5 f 1 

y = a^'^ cosh~^ < J — z } (27) 

^ 1874161 IV 74 J ^ ^ 

Let us consider the case of n = |. As this takes place Og^ have the 
following values 

= 1, a[,^^ = 2, = 3, a[,^^ = -3 (28) 

We see that elliptic solution may be obtained if we consider solutions 
corresponding to combination of a^^^ = 3 and a^^^ or Cq^^ Oq^^ and 

Consider the solution correspond to coefficients a^^^ = 3 and a^f^ . Accord- 
ingly to method from papers [T],|2] we construct the formal Laurent expansion 
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of solution for Eq. ([8]) at n = | corresponding a^f^^ = —3 



3 a f S 2lla^ \ , 
w = —z - —— + z^- 



z 34 V210 242760 



(29) 



a (679 + 3468 5) 5 9 

— Z -\- . . . -\- Oj-\r\ Z -|-... 

24761520 ^ 10 ^ 

Note that coefficient aio in series fl29|) is arbitrary constant. 
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Following [T1I2] we look for the elliptic solution of Eq. in the form 

Comparing the expansion (129|) with the Laurent expansion for expression 
fl5U]) we have 

a , 225 

^0 = 0, A= , 5 = 0, 6 = , 

" ' 102' ' 1156 ' 

708451848720' 132651' 867' 



aiQ 



708451848720 



We see that there is constraint on the coefficient aio- Otherwise solution 
in the form f[T7|) does not exist. 

In this case the Weierstrass elliptic function is degenerated to trigono- 
metric function as well 



p(z - zo, — , ) = — ( 3 tan^ I J — iz -zq) \ - 2] (32) 

°'867' 13265r 102 1 11/34^ °^ ' ' ^ ' 



After substituting formulae fl32|) into fl30l) we have solution of ([8]) in the 
form 

Using the Backlund transformations ([7j) we obtain exact solution of the 
Eq.([6]) in the form 



y 



90 \ 
289/ 



3/4 



289 / V 34 



a3/2cosh-3<|,/-^4 (34) 




In the case of n = 7/3 and solution corresponding to combination to three 
asymptotic expansions ag^^ ag^^ Oq'^'* exactly the same solution as (133|) . 

So we see that Eq. (|6]) has the elliptic solutions only at n = 2, n = 3 and 
n = 5. This completes the prof. □ 

Remark 1. Let us consider the case of a = separately. Comparing the 
expansion (12 9 p at cr = with the Laurent expansion for function (130 p in the 
case of n = I we have 

^3 ^ 

A = S = /lo = 0, aio = , 92 = , 93 = (35) 

" ' " 2167074000' ^ 63' ^ 
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Using this relation we obtain the eUiptic solution of Eq. in the form 

3p, (^,-^,0) 

w = 7 ^ — :r 36 

The case of cr = means that Cq = a. So we can see that elliptic solutions 
exist only at special value of Cq. 

Using the Backlund transformations ([7]) from (15^ we obtain the elliptic 
solution of Eq. (jS]) in the form 



p(^,-^,0) 



In the cases n = ^ and n 



3/2 
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(37) 



g at cr = elliptic solutions are degenerated 
to the rational solutions under condition 6 = 0. 

Remark 2 The simple periodic solution of Eq. has the following form 

47r 



w 



cot — - 



- 1) T V T 
where T and a have values 



(3^ 



-2v^(n + l)7r 

(n- l)v^ ' 
5{n^ + 2n + 5) 



0^1,2 



2v^(n + 1) 



3,4 



0"3,4 



2i/2v^(r2 + l)7r 

±— ^ — , 

(n-l)v^ 

^(n^ + 2r2 + 5) 



(39) 



2V^(Ti + 1) 



Using formulae (138|) . relations ( 139|) and taking Backlund transformations ([7]) 
into account we get 

/ \^ 

(3^2 + lOn + 3)5 



1/1,2 



v 



+ 1) 



cos 



V=F(«+l)v^ 

Expression fHUl) is the simple periodic solution of Eq 



(40) 



at any n 7^ ±1. 



3 Exact solutions of Eq. (H) in the case of 
n = 5 and m 7^ 1. 

Suppose that solutions of Eq. 02]) have the asymptotic expansions in the 
form of the Laurent series in a neighborhood of the pole z = Zq 

oo 

y{z) = Y,(^k{z-z,f-P, p>0 (41) 

A:=0 
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Eq.(|n]) is autonomous and without loss of the generahty we can assume 
^0 = 0. 

In the case of n = 5 Eq. ([3j) has the four different formal Laurent expan- 
sions in a neighborhood pole 2; = of the first order {p = 1). 

Without loss of generality we use the following values of parameters 7 
and f3 

/3 = -l, 7 = 24 (42) 

Corresponding values of the coefficients Oq are the following 

a« = 1, = -1, = I, a^o^ = -%. (43) 

The necessary conditions for existence of the elliptic solutions take the 
form 

af + af = and a'i^ + a^^^ = (44) 

We see that the coefficients a'^^ and Oq'^-' correspond to the elliptic solu- 
tions. 

The Fuchs indices corresponding of each expansions are following 
Ji = -1, J2,3 = ^(5±v^) (45) 

So we see that coefficients of expansion fHTj) for the solution of Eq. ([3]) are 
unique. 

We use the method from the references |Ti,[2] for constructing the elliptic 
solutions Eq. The algorithm of this method is the following. 

1) Construct the formal Laurent expansion of solution for Eq.dH]); 

2) Take the general form of the possible elliptic for Eq.(|3]) and find the 
formal Laurent expansion for this form of solution; 

3) Compare the formal Laurent expansion of solution for Eq.([3]) hat was 
found on the first step with the formal Laurent expansion of the possible 
eUiptic solution that was found on the second step; 

4) Solve the system of the algebraic equations obtained on the third step 
and find the parameters of the Eq. fH6|) and parameter of the possible elliptic 
solution. 

Remark 1. The formal Laurent expansion for general form of the possi- 
ble elliptic solution can be found by using textbook [IB] or by using Maple. 

Remark 2. The simple periodic solution can be found by using analogues 
algorithm. See details in references pljl2]- 

Consider the exact solutions of Eq. ([3]) in the case m = 2. 
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Substituting m = 2 into Eq. in the case of n = 5 we have 

yl - {y-z y... - \ yl) + {y' + ^y' + c, = ^ (46) 

The formal Laurent expansion of the solution for Eg. (1461) corresponding 
to Oq^^ = 1 is the following 

^ z \ m ) 120 1800 5760 

The possible exact solution of Eq. (j3H]) takes the form 

y= 2{p-A) 

where p is Weierstrass elliptic function. 

Comparing the expansion Wl\ with the Laurent expansion of expression 
IHl) we have 

/io = 0, A = 0, B = ^, Co = ±v^, 

92 = 0, Os = , C2 = 

^ ' 3600' 1800 

Using these relations we obtain the elliptic solution of Eq. (H6|) in the form 



P {z, 0, ^ ^^^) + 
2p(2;,0, 3goo, 



^ = -L^-1__3600Z^ (50) 



Let us look for the exact solutions of Eq. in the case of m = 3. 
Substituting m = 3 into Eq. ([3]) in the case of n = 5 we have 

yl - {y. y... - \ yl) + {y' + Ay' + c, = q (51) 

The formal Laurent expansion of the solution for Eq. (l46i) corresponding 
to Oq^'' = 1 is the following 

1 ( 2{a-Cl)~5 \ _ {2{Cl-a) + 5) [A{Cl - a) - 5) 
^ " z V 120 ) 14400 ^ ■ ■ ■ 

(52) 
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Comparing the expansion (15^ with the Laurent expansion for the expression 
(jH]) we have 

120 ' ' ^ 240 

_ (13 + 13 a) (-2 g + 2 - 6f .^gx 
~ 432000 

^ (2(Cg -a) -6) (8(Cg - g) + ^) (18(Cg - a) + ^) 
^ 96000 

Using these relations we obtain the elhptic solutions of Eq. (15T]) in the form 

y = p'jz- zo,g2,g3) ^^^^ 

where 92,93 are determined by the relations (153|) 

Let us look for the periodic solutions of Eq. (!5T|) . In accordance with the 
method from the references [Ul^] we search for the periodic solutions in the 
form 

y = -^fl- V \ rr ^ + v^cot(v^2;) + (55) 

A - \'L coi{\/ Lz) 

Comparing the expansion f l52|) with the Laurent expansion for the possible 
solution f l55|) we have 

5 = 3(a-Co^), L = ^^^ A = V^, ^^^^ 
ho = —^/—L, C2 = —L^ 

Let us consider the exact solutions of Eq. in the case of m = 4. 
Substituting m = 4 into Eq. in the case = 5 we have 

yl - (y^ I/... - 1 yl) + ^ + 4 / + = o (57) 

The formal Laurent expansion of solution for Eq. fl46p corresponding to coef- 
ficient Cq^'' = 1 is the following 



1 6 /120(a-C2) + 52\ 53 

y = h ^ — z z^ + ... (5^ 

^ z 120 V 7200 / 432000 ^ 
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Comparing the expansion (15^ with the Laurent expansion for the expression 
we have 

A = -,T,r:rT, B = ho - 



43200' 216000' 120' 

~ ~ 31 104000' ~ 1259712000000' ^ 



^ 186624000000' ° 12 V 5 
Using these relations we obtain the elhptic solutions of Eq. fH6|) in the form 

y ^ - ZQ, 92,93) + 21I0OO _ ^ ^gQ^J 

where 92,93 are determined by relations fl59|) . 



4 Exact solutions of Eq. (Oj) at m = 2 in the 
case n = 3. 

Let us consider the exact solutions of Eq. in the case of n = 3 and m = 2. 
The case of n = 3 and m = 1 was considered in [2]. From Eq. we have 

y.^ + I/... - I yl) + ly' + ly' + C2 = o. (6i) 

Without loss of the generality we take the following values of /3 and 7 
/3 = -1, 7 = 120. (62) 

We use the algorithm again from the reverences [Tl[2] that we have de- 
scribed in the previous section. 

We assume that the solution of Eq. fl6T]) can be presented in the form of 
the Laurent series f HT]) . Without loss of the generahty we set zq = 0. 

Eq. fl6T]) admits two different Laurent expansions in the neighborhood of 
the pole z = 



_ 1 , 6(Cg -a) -5 (-^ + 6Cg-6a)(^ + 6Cg-6a) , , 

360 129600 ^+...(63j 



1 6{a-Ci)-S (-S + 6Ci-6a){6 + 6C'o-6a) , 
z^ 360 129600 "' ^ ' 
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We see that aL^j and a^] are zeros. So the necessary condition for the 
existence of the elhptic solution is hold. 

We use the expression for the solution of Eq. fl61l) in the following form 



y = 2^(2,^2,^3) - 



p' {2,02,93) +B 
p{z,92,93) - A 



+ ho + A 



(65) 



Comparing the expansion (163|) with the Laurent expansions to (l65l) we 
have 



60 



Co 



540 
303750 



93 



{CI - a) 
81000 



(66) 



After substituting this relations in formulae dHS]) we get following solution 
of Eq.dSU 



y = 2p U, 



{Cl-af {CI -ay 



540 



81000 



1 

'4 



P 1 ^5 540 ' 81000 



,2 /„o \3\ ^2 
'P I ^' 540 ' 81000 j ~' 60~ 



+ 



a-Cl 
60 



(67) 



We should note that this solution can be obtained from solution (4.19) in 
the reference [2] assuming Co = 0. 



5 Conclusion 

In this paper we have studied all possible solutions of the generalized Brether- 
ton equation. The meromorphic exact solutions of this equation are classified 
in the cases n = 5atm = 2, m = 3 and m = 4. The case m = 1 and ar- 
bitrary n also were considered. We have shown that the elliptic solutions of 
the generalized Bretherton equation exist only in the cases n = 2, n = 3 and 
= 5 at 0" 7^ and at 6 ^ 0. We have obtained that there is the elliptic 
solution in the case n = | at Cq = a. 
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